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Abstract 

Decays of the type B -» D^D^ can be used to provide a measurement of the parameter sin2/3 of 
the Unitarity Triangle that is complementary to that derived from the mode B° —* J/tp K® . Here we 
report a measurement of the branching fraction and a study of the CP parity content for the decay 
B° — > D* + D*~ with the BABAR detector. With data corresponding to an integrated luminosity of 
20.7 fb" 1 collected at the T(45) resonance during 1999-2000, we determine the branching fraction 
to be B(B° -> D* + D*-) = (8.0 ± l.6(stat) ± 1.2(syst)) x 10~ 4 . The measured fraction of the 
component with odd CP parity is 0.22 ± 0.18(stat) ± 0.03(syst). Observation of a significant 
number of candidates in the decay modes B° — ► D* + D~ and B + — > D*+D*° [ s reported. All 
results presented in this note are preliminary. 
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1 Introduction 



One of the most important goals of the BABAR experiment is to precisely measure the angles of 
the Unitarity Triangle. While the decay B° — > J/ip can be used to measure sin2/3, the Standard 
Model predicts that the time-dependent CP-violating asymmetries in the decays B° — > 
can also be used to measure the same quantity |IJ. An independent measurement of sin2/3 in 
B° —* Z)(*)+I)(*)~ modes is especially important since several typical extensions to the Standard 
Model can lead to different assymmetries between this quantity in B° —* Z)(*)+Z)(*) _ events and 
charmonium events ||). Measurements of sin2/3 in B° — > £)(*)+£)(*) - would thus provide stringent 
tests of CP-violation in the Standard Model and have significant potential to indicate deviations 
from it. However, the vector- vector decay B° — > D* + D*~ is not a pure CP eigenstate and a 
sizeable dilution of the measured asymmetry can be produced by a non-negligible P-wave CP-odd 
component. The dilution can, in principle, be completely removed by a time-dependent angular 
analysis of the decay products |J. 

The rate for the Cabibbo-suppressed decays B — > £)(*).£)(*) can be estimated from the measured 
rate of the Cabibbo-favored decays B — > D^D^: 

B{B -> DWflW) « ( ] tan 2 6 C ■ B(B -» D^D^), (1) 

where 0^ is the Cabibbo angle, and /^(.j and /„(*) are decay constants. From this it follows 
_ u s 

that the B — > £)(*)£)(*) branching fractions are of the order of 10 -3 . Previous measurements of 
branching fractions and upper limits for these modes are summarized in Table |l[ 

In section |2| we describe briefly the BABAR detector and the dataset. Section ^ describes the 
measurement of the branching fraction B° — > D* + D*~ and section Q the corresponding angular 
analysis to extract the CP-even component of the decay. Finally, in section |5| we discuss the 
observation of the decays B° -> D* + D~ and B+ -» D*+D*°. Section |6| summarizes our results. 

Table 1: Summary of branching fraction and upper limit measurements performed by the CLEO 
Q and ALEPH || experiments. Upper limits are quoted at the 90% confidence level. 





Decay 


Branching Fraction (x 10 4 ) 


B u 


-> D*+D*- I 


| 9.9l^(stot)±1.2(syst) 


B° 


-> £)*+£)- [ 


| <6.3 


B+ 


-»■ D*+D*° f 


| < 110. 



2 The BABAR detector and dataset 

The data used in this analysis were collected with the BABAR detector || at the PEP-II storage 
ring Q located at the Stanford Linear Accelerator Center. This data sample represents an inte- 
grated luminosity of 23.3 fb — , with 20.7 fb _1 collected on the T(AS) resonance. The total number 
of BB pairs produced in this sample was N B g = (22.7 ± 0.36) x 10 6 . 

Charged particles are detected and their momenta measured with the combination of a 40-layer 
drift chamber (DCH) and a five-layer silicon vertex tracker (SVT) embedded in a 1.5 T solenoidal 
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magnetic field. Photons are detected by a Csl electromagnetic calorimeter (EMC) that provides 
excellent angular and energy resolutions with a high efficiency for energies above 20MeV. Charged 
particle identification is provided by the specific ionization loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering the barrel region of 
the detector. 



3 Measurement of the B° — > D* + D* branching Fraction 

B° mesons are exclusively reconstructed by combining two charged D* candidates reconstructed in 
a number of D* and D decay modes. Events are pre-selected by requiring that there be three or 
more charged tracks and that the normalized second Fox- Wolfram moment || of the event be less 
than 0.6. We also require that the cosine of the angle between the reconstructed B direction and 
the thrust axis of the rest of the event be less than 0.9. 

Charged kaon candidates are required to be inconsistent with the pion hypothesis, as inferred 
from the Cherenkov angle measured by the DIRC and the specific ionization measured by the 
SVT and DCH. No particle identification requirements are made for the kaon from the decay 

D° K'TT + . 

K® — » 7r + 7r~ candidates are required to have an invariant mass within 25 MeV/c 2 of the nominal 
Kg mass ||. The opening angle between the flight direction and the momentum vector of the K® 
candidate is required to be less than 200 mrad, and the transverse flight distance from the primary 
event vertex must be greater than 2 mm. 

Neutral pion candidates are formed from pairs of photons in the EMC with energy above 30 MeV, 
an invariant mass within 20 MeV/ c 2 of the nominal 7r° mass, and a summed energy greater than 
200 MeV. A mass-constraint fit is then applied to these 7r° candidates. The ir° from D* + —* D + ir° 
decays ("soft" n ), however, is required to have an invariant mass within 35 MeV/ c 2 of the nominal 
7T° mass and momentum in the T(4S) frame of 70 < p* < 450MeV/c, with no requirement on the 
summed photon energies. 

The decay modes of the D° and D + used in this analysis were selected by an optimization of 
S 2 /(S + B) based on Monte Carlo simulations, where S is the expected number of signal events and 
B is the expected number of background events. The D° and D + modes used and their branching 
fractions are summarized in Table ||. D° (D + ) meson candidates are required to have an invariant 
mass within 20 MeV/ c 2 of the nominal D° (D+) mass. 

The D* + mesons are reconstructed in their decays D* + — ► D°ir + and D* + — > D + tt°. We include 
for this analysis the decay combinations D* + D*~decaying to (D°7r + , D°ir~) or (D°tt + , D~n°), but 
not (D + 7T°,-D - 7r ) due to the smaller branching fraction and larger expected backgrounds. The 
branching fractions for these modes are summarized in Table ||. D° and D + candidates are subjected 
to a mass-constraint fit and then combined with soft pion candidates. A vertex fit is performed 
that includes the position of the beam spot to improve the angular resolution of the soft pion. 

To select B° candidates with well reconstructed D* and D mesons, we construct a x 2 that 
includes all measured D* and D masses: 



2 _ / "'/> - i"Pri,G \ 2 , ( m D- m Dp D g 

\.!l-!.vs 1 I ) r I ~ 

/ Amn* — Amn* \ 2 / Am-^* — Amn- 

' ' U PDG III O U PDG 



&Am D , J \ CTA 
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Table 2: D° and D + decay modes and branching fractions 0. The branching fraction for 
Kg — > 7r + 7r _ is included for modes containing a ifg. 



Decay Mode 




Branching Fraction (%) 


D" -► K--K+ 

D° -► K-1T+1T 

D° -> K^vr+Tr+vr" 




3.83 ± 0.09 
13.9 ±0.9 
7.49 ± 0.31 
1.85 ±0.14 


Total D u Branching 


Fraction 


27.1 


Decay Mode 




Branching Fraction (%) 






9.0 ±0.6 
0.99 ± 0.09 
0.87 ±0.07 


Total D + Branching 


; Fraction 


10.9 



Table 3: D* and D*° decay modes and branching fractions |J. D*° is used for the B + — > D* + D* { 
analysis described in Section |5[ 



Particle 


Decay Mode 


Branching Fraction (%) 


D *+ 


D*+ -» 


J D U 7T+ 


67.7 ±0.5 




D*+ -» 


D+vr 


30.7 ±0.5 


Total Visible Z?* + Branching Fraction 






98.4 




L>* u -» 


L> u vr u 


61.9 ±2.9 




D*° - 




38.1 ±2.9 


Total D*° Branching Fraction 






100.0 



where the subscript PDG refers to the nominal value, and Am/)* is the D* — D mass difference. 
For a rriD we use values computed for each D candidate, while for u A.m D * we use fixed values of 
0.83MeV/c 2 for D* + -» D°tt+ and 1.18MeV/c 2 for D* + -» D + tt°. A requirement that xltass < 20 
is applied to all B° candidates. In events with more than one B° candidate, we choose the candidate 
with the lowest value of XMass- 

A B meson candidate is characterized by two kinematic variables. We use the energy-substituted 
mass, mES) defined as 

m E s = \/ E* B 2 eam - p* B 2 
and the difference of the B candidate's energy from the beam energy, AE , 

AE = E* B - E* Beam 

where E* B {p* B ) are the energy (momentum) of the B candidate in the center-of-mass frame and 
E* Beam is one-half of the center-of-mass energy. The signal region in the AE vs. ities plane is defined 
to be \AE\ < 25MeV and 5.273 < m^s < 5.285 GeV/c 2 . The width of this region corresponds to 
approximately ±2.5cr in both AE and twes- 
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These values on Xm„ ss1 "1es> and AE were chosen based on an optimization of S 2 /(S + B), 
where S is the expected number of signal events and B is the expected number of background 
events. The optimization process was done entirely with samples of signal and generic BB and cc 
Monte Carlo where the background distribution is taken from a sideband region, defined as 

\AE\ < 200 MeV 

5.20GeV/c 2 < m ES < 5.26GeV/c 2 

and 

50 MeV < |A£7| < 200 MeV 

5.26 < m ES < 5.29GeV/c 2 

These values were chosen based on a maximization of S 2 /{S+B) with a tendency towards looser cut 
values to reduce any possible systematic error incurred due to the differences in the reconstructed 
mass resolutions between data and Monte Carlo. 

To determine the number of signal events in the signal region, we must estimate the expected 
contribution from background. This is done by scaling the number of events seen in the data 
sideband (defined above) with a scaling factor which gives a measure of the relative areas of the 
signal region to the sideband region. We parameterize the shape of the background in the AE vs. 
tties plane as the product of an ARGUS function [1C] in m^s and a first-order polynomial in AE . 



Based on this parameterization we estimate that the ratio of the number of background events in 
the signal region to the number in the sideband region is (1.72 db 0.10) x 10~ 2 . The uncertainty is 
derived from the observed variation of this ratio under alternative assumptions for the background 
shape in tbes and AE . 

Figure |l| shows the events in the AE vs. toes plane after all selection criteria have been applied. 
The small box in the figure indicates the signal region defined above, and the sideband is the entire 
plane excluding the region bounded by the larger box outside the signal region. There are a total of 
38 events located in the signal region, with 363 events in the sideband region. The latter, together 
with the effective ratio of areas of the signal region to the sideband region, implies an expected 
number of background events in the signal region of 6.24 db 0.33(stat) db 0.36(syst). The quoted 
systematic uncertainty comes from the background shape variation discussed previously. Figure |2| 
shows a projection of the data on the mES axis after requiring \ AE\ < 25 MeV. 

We use a detailed Monte Carlo simulation of the BABAR detector to determine the efficiency for 
reconstructing the signal. This, together with the total number of BB pairs produced during data 
collection, allows us to determine a preliminary branching fraction for B — > D* + D*~ to be 

B{B° -> D* + D*-) = (8.0 ± l.6(stat) ± 1.2(syst)) x 10" 4 

The dominant systematic uncertainty in this measurement comes from our level of understand- 
ing of the charged particle tracking efficiency (9.4%). The high charged particle multiplicity in this 
decay mode makes this measurement particularly sensitive to tracking efficiency. Systematic errors 
were assigned on a per track basis for tt, K and slow it, and were added linearly due to large cor- 
relations. The imprecisely known partial-wave content of the B° — > D* + D*~ final state is another 
source of systematic biases. Monte Carlo events in each of the two extremes of transversity ampli- 
tudes (A // ,V2A ,A±) = (l.,0.,0.) and (0.,1.,0.) were generated and reconstructed [[H]] . Although 
both mixtures correspond to Rt = 0, the resulting pt distributions of the slow pion represent the 
two extreme cases of possible pt distributions. The change in the reconstruction efficiency of these 
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Figure 1: Distribution of B° — ► D* + D* events in the AE vs. toes plane. The small box indicates 
the signal region, while the sideband region is everything outside the larger box. 
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Figure 2: Distribution of B° — > D* + D* events in toes plane with a cut of \ AE\ < 25 MeV applied. 
The curve represents a fit to the distribution of the sum of a Gaussian to model the signal and an 



ARGUS function [1C] to model the background shape. 



final angular states is quoted as systematic error (6.6%). Other significant systematic biases arise 
due to the uncertainties on the D* + , D° and D + branching fractions (5.6%) and the differences 
in mass resolutions between Monte Carlo and data (4.1%). Possible contributions from peaking 
backgrounds was found to be negligable.The total systematic uncertainty from all sources is 14.5%. 
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4 Determination of the CP content of B° — ► D* + D* 

The fraction of the CP-odd component, Rt, of the decay B° — ► D* + D*~ can be determined from 
the angular distribution in the transversity basis ||: 

^ T ^V = 7( 1 -^) sin2 ^ + ^ cos2 ^ ( 2 ) 
T d cos 0tr 4 2 

Here T is the decay rate and 9t r is the polar angle defined as the angle between the normal to the 
D*~ decay plane and the tt + line of flight in the D* + rest frame. 

We perform an unbinned maximum likelihood fit to the 38 events in the signal region described 
in the previous section. The fit takes into account the presence of background, whose properties 
are derived from the sideband sample, and the angular resolution ag tr estimated from the simulated 
data samples. We define the likelihood function to be 

C= Y[Ci= J] [pxT{Otr,i,oe,uR S t g )+ 

8=1, n i=l,n ^3) 

(l-p)xF(e tr4 ,a e4 ,R? 9 ) 

where n is the number of selected events and the contribution to the total likelihood from the i-th 
event, Ci, is defined in terms of the purity p of the sample and the probability density functions 
^ : (&tr,i,o'e,i, Rt) for the signal and background: Rf 9 and Rf. 9 are the parameters describing the 
shapes of the signal and background angular distributions, respectively, and 9t r ,i is the measured 
transversity angle in event i. The probability density functions T are obtained from the convolution 
of the angular distribution (Eq. |2|) with Gaussian resolution functions describing the measurement 
errors a$^. We note that Eq. |2| is the differential decay rate T integrated over the full ranges of 
the other two decay angles in the transversity basis; this assumes a flat acceptance and needs to 
be corrected for in the final determination of Rt- An empirical study of the events in the sideband 
region showed that the background angular distribution can be described as constant in cos 6t r , as 
expected. The parameterization in Eq. |||can describe deviations from this behaviour and was used 
to estimate the background contribution to the systematic uncertainty on Rt. 

The fit procedure was tested on several samples of B° — > D* + D*~ Monte Carlo events generated 
with different Rt values and different sets of decay amplitudes in order to determine the possible 
bias induced by angular acceptance effects correlated with the inefficiency in detecting soft pions 
from D* decays below a threshold in transverse momentum of about 70 MeV/ c. The fitted Rt values 
were fully consistent with the generated values in the limit of negligible soft pion inefficiency, but 
could be biased, depending on the decay amplitudes, when the pion detection threshold was taken 
into account. A correction and a systematic error were estimated to describe these effects (described 
below) . 

The value of R hkg was evaluated by fitting the 363 events in the sideband region and setting 
p = in Eq. ||. The result of this fit was R bk9 = 0.29 ± 0.04, compatible with the value expected 
for a flat distribution (Rt = 1/3). 

To determine R s t 19 , we fit the 38 events in the signal region with R h t k9 fixed to 0.29 and p fixed 
at 83.6%. The result of the fit to the signal region, without correction for angular acceptance bias, 
is R s t 19 = 0.25 db 0.18(stat), and is shown in Figure |3[ The fit was repeated with both R s t 19 and 
R h t k9 floating, giving the same central values with a rather small correlation (—0.04) between the 
parameters. 
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Figure 3: The cos 6t r distribution from the unbinned ML fit, superimposed on the histogram of 
the B° — > D* + D*~ candidates in the signal region. The dotted and dashed lines represent the 
fitted CP components for the signal and background contributions. The lower two dotted curves 
represent the two background fitted components, proportional to sin 2 9t r and cos 2 6t r respectively; 
their sum (dashed curve) is almost constant at about 0.3 events/bin. Similarly, the two upper 
dashed curves represent the CP components of the signal. The solid line represents the cos 6t r 
distribution from the unbinned ML fit for the selected events. 



The dominant systematic uncertainty on the measured value of Rt comes from the acceptance 
bias due to the transverse momentum threshold in the detection of soft pions. This bias ranges 
between —0.048 and +0.004, and depends on the yet unknown decay amplitudes. The central value 
of this interval was taken as a correction to the fitted Rf 9 , while its half width was assumed as a 
contribution to the systematic error. Other contributions come from the imperfect knowledge of the 
finite resolution on the measured value of 9t r (0.006), of the angular distribution of the background 
(0.008) and of the purity of the signal sample (0.0003). All uncertainties were evaluated with Monte 
Carlo simulations. The total systematic uncertainty on Rt was determined to be 0.031, giving the 
final corrected result: 

R t = 0.22 ± 0.18(stat) ± 0.03(syst) 

The likelihood function was also used to estimate approximate confidence intervals, taking into 
account the allowed physical region. Fig. |] shows the dependence of L^taiRt) = — 21n(£(i? t )/£(fii)) 
on Rt, where Rt is the value of Rt that maximizes C(Rt). The lines corresponding to the 68%, 
95%, 99% confidence level limits on R t are also plotted. We can exclude values of Rt greater than 
0.63 at the 95% confidence level. 
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Figure 4: Likelihood ratio Ldata(Rt) = — 21n(£(i2t)/£(^)) for the data, where Rt is the value of 
Rt that maximizes C(Rt), and 68%, 95%, 99% confidence level limits on R t . The systematic error 
on Rt is indicated by the dotted and dashed lines. 



5 Observation of the decays B° -> D* + D and B + -> D* + D*° 

The decays B° — > D* + D~ and — > D*+D*° are studied following a method largely similar to 
that described in Sec. ^. Here, only those aspects of the analyses that differ significantly from that 
of the B° -> D*+D*- analysis are discussed in some detail. 

For B° — > D* + D~, B° mesons are exclusively reconstructed by combining a Z>* ± and a D T 
candidate that are reconstructed in a number of D*^ and decay modes. For B + — ► 
the exclusive reconstruction combines a Z}* 1 * 1 and a £)*°. The kaon flavor of the D*° is checked 
to make sure that a D* + is paired only with a D*° and a D*~ is only paired with a The 
selection of and D* 1 * 1 candidates, and the K® and 7r° candidates that are used to compose them, 
is identical to that described for the B° — ► D* + D*~ analysis. 

The decay modes of the D and D* used in these analyses are selected by an optimization of 
S 2 /(S + B) based on Monte Carlo simulations. D*^ mesons are reconstructed in their decays 
D* + — > Z?°7r + and D* + — ► L> + 7r°, and Z?*° mesons are reconstructed in their decays — > D°tt° 
and -D* — > D°7. Modes used and their branching fractions are summarized in Tables ^ and [|. As 
in the B° — > D* + D*~ analysis, we construct XMass variables that include all measured Z)* ± , D*°, 
and D masses. For B + — > D* + D*°, XMass contains 4 terms: 
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For B — > D* + D , XMass contains 3 terms: 

2 _ / m D - m DpDG \ 2 / m Dn , - m DpDG 

AMass — \ I "t" 

V <>Vn D / y %d d , 

The major difference between these analyses and the B° —* D* + D*~ analysis is that the XMass 
cut values for these analyses are set individually for each submode instead of having a global XMass 
value for all submodes, to better take into account the fact that the amount of background is quite 
different for each of the different submodes in these analyses. The procedure used to optimize the 
XMass value chooses the individual values of XMass that together maximize the global S 2 / (S + B). 
Submodes which have an optimal value of less than 2 times the number of degrees of freedom of the 
XMass (8 f° r B + — ► D* + D*°, 6 for B° — > D* + D~) are rejected on the grounds that the resulting 
sensitivity is poor, and that the tightness of the XMass vame makes such modes more susceptible 
to systematic errors in measured yield. 

These values were tuned with samples of signal and generic BB and cc Monte Carlo where 
the background distribution is taken from a sideband region. For B + — > D* + D*°, the sideband 
region is the same as that for B° — > D* + D*~ . For B° — > D* + D~ , the region where decays (such as 
B° — > D* + D*~) can feed down into the AE -tbes plane must be eliminated from the sideband. As 
B° — > D* + D*~ contains a reconstructed D*^ and D, it is only separated from B° — > D* + D~ events 
due to the missing energy of the slow pion from the second D*^. This missing energy manifests 
itelf as a negative shift in AE . These events accumulate in the area below the B° -» £>*+.D- 
signal region, and in order to remove them, the region defined by: 

AE < -50MeV 

m ES > 5.26 GeV/c 2 

is removed. The toes and AE distributions for events reconstructed in the channels B° — ► D* + D~ 
and B + — > D* + D*° are shown in Fig. |5| and in Fig. |6|. 

For the -> £>*+£>" channel we reconstruct a total of 31 events, of which 10.5 ± 1.7(stot) are 
background. The probability that the visible signal is a statistical fluctuation of the background 
is 9.7 x 10~ 7 (> 4.3a). As pointed out in the introduction, this channel is a CP conjugate state 
that can be used for sin2/3 measurements. For B + — > D* + D*° we reconstruct a total of 39 events, 
of which 20.3 ± 0.5(stat) are background. The probability that the visible signal is a statistical 
fluctuation of the background is 2.9 x 10~ 6 (> 4.1a). This channel is useful for calibrations and 
tests of the CP fitting procedures. 



6 Summary 

Using data collected by the BABAR experiment during 1999-2000, we have observed a signal of 
31.8 ± 6.2(stat) ±0.4(syst) events in the decay B° — > D* + D*~ . Our measurement of the branching 
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Figure 5: Top Left: w-es projection of the B — > D* + D~ event population, in the Ai? signal 
band (—0.025 < AE < 0.025 GeV). The crosses are the data. The dashed line represents the 
extrapolation, to the AE signal band, of a two-dimensional background fit in the AE -m^s 
sidebands. The solid line is the sum of this background extrapolation and of a fitted, Gaussian- 
shaped signal centered on the B mass. Top Right: AE projection of the B° — ► D* + D~ event 
population, in the mES signal band (5.273 < uies < 5.285 GeV/ c 2 ). The crosses are the data. 
The dashed line represents the extrapolation, to the tties signal band, of the above-mentioned 
background fit. The solid line is the sum of this background extrapolation, and of a gaussian- 
shaped signal centered on AE = 0. In the case of both the upper plots, the points in the feed-down 
region (AE < —0.05 GeV) are excluded from the fits. Bottom: Two-dimensional distribution of 
the B° — > D* + D~ events in the AE vs. tues plane. The small ellipse indicates the signal region, 
while the sideband region is everything that is outside the box that surrounds the signal region and 
also outside the (feed-down) box below the signal region. 



17 



> 

u 
O 
in 

O 

o 
d 



> 

UJ 



PRELIMINARY 




> 
o 

o 
© 



> 



5.24 5.26 5.28 5.3 

Energy-Substituted Mass (GeV) 




BaBar 



0.05 



0.1 0.15 0.2 
delta E (GeV) 



0.2 



BaBar 



o.i 



> 
a 



-0.1 



-0.2 



5.2 




5.26 5.28 
m ES (GeV/c 2 ) 



Figure 6: Top Left: toes projection of the i3 + — > D* + D*° event population, in the AE signal 
band (—0.025 < AE < 0.025 GeV). The crosses are the data. The dashed line represents the 
extrapolation, to the AE signal band, of a two-dimensional background fit in the AE -toes 
sidebands. The solid line is the sum of this background extrapolation and of a fitted, Gaussian- 
shaped signal centered on the B mass. Top Right: AE projection of the B + -> D*+D*° event 
population, in the toes signal band (5.273 < toes < 5.285 GeV/ c 2 ). The crosses are the data. 
The dashed line represents the extrapolation, to the toes signal band, of the above-mentioned 
background fit. The solid line is the sum of this background extrapolation, and of a Gaussian- 
shaped signal centered on AE = 0. Bottom: Two-dimensional distribution of the B + — > D* + D*° 
events in the AE vs. toes plane. The small ellipse indicates the signal region, while the sideband 
region is everything that is outside the box that surrounds the signal region. 
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ratio is 

B{B° -► D*+D*-) = (8.0 ± l.Q(stat) ± 1.2(syst)) x 10" 4 

From the transversity angular distribution of these events, the measured fraction of the component 
with odd CP parity is 

R t = 0.22 ± 0.18(stat) ± 0.03(syst) 
Finally, signals are also observed in the decay modes B° — ► D* + D~ and S + — > 
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